arXiv:1505.03879vl [gr-qc] 14 May 2015 


Superfluidity as a Tool for Constraining the Energy—Momentum Relation 


A. Camacho* and S. Guti&rez^ 

Departamento de Fisica, Universidad Autonoma Metropolitana-Iztapalapa 
Apartado Postal 55-534, C.P. 09340, Mexico, D.F., Mexico. 

(Dated: May 18, 2015) 

The present work explores the possibilities that superfluidity could offer in the context of quantum 
gravity phenomenology, at least in the realm of deformed dispersion relations. The experimental pro¬ 
posal involves a Bose-condensed sodium gas trapped by an isotropic harmonic oscillator potential. 
A deformed dispersion relation for the particles of the system is considered and the consequences 
of this assumption upon the possible region of superfluidity of this system is analyzed. It will be 
shown that in this sense the effects of quantum gravity could be detected resorting to experiments 
of superfluidity in Bose-Einstein condensates. Finally, using the current experimental results in this 
direction an upper bound for the corresponding phenomenological parameters will be also obtained. 
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I. INTRODUCTION 

Quantum gravity phenomenology [1-3] emerges as an 
answer of the community to the mathematical and phys¬ 
ical difficulties plaguing all theoretical models behind a 
quantum theory of gravity [4, 5]. These efforts embody 
new physical effects, for instance, deformed versions of 
the dispersion relation, deviations from the 1 /r-potential 
and violations of the equivalence principle. At this point 
it is noteworthy to comment that these cases do not ex¬ 
haust the extant possibilities. 

In the quest for a solution to this long-standing puz¬ 
zle in modern physics we have some efforts which entail, 
unavoidably, the breakdown of Lorentz symmetry [6-8]. 
Lorentz symmetry is a the bedrock of modern physics 
and, therefore, it has been subjected to some of the high¬ 
est precision tests in Physics [9-11]. The current exper¬ 
imental results show no evidence of a violation of this 
symmetry, nevertheless this last fact does not discard it 
and, in consequence, further work is required. At this 
point it has to be clearly stated that the phrase violation 
of Lorentz symmetry has several meanings, i.e., it embod¬ 
ies several characteristics. For instance, Local Lorentz 
Invariance, or Local Position Invariance [12]. For us it 
will mean a modification of the dispersion relation. Let 
us state this phrase in a mathematical language. As men¬ 
tioned above several quantum-gravity models predict a 
modified dispersion relation [6-8], the one can be charac¬ 
terized, phenomenologically, through corrections hinging 
upon Planck’s energy, Ep 





( 1 ) 


Here a is a coefficient, whose precise value depends 
upon the considered quantum-gravity model, while n, 
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the lowest power in Planck’s length leading to a non¬ 
vanishing contribution, is also model dependent. 

The quest in this direction has already considered in¬ 
terferometry as a tool [13-17]. Clearly, all these previous 
works in quantum gravity phenomenology involve differ¬ 
ent fields of physics [18], i.e., the search includes now 
many areas of modern physics. One of these topics is 
condensed matter physics. Indeed, the use of cold atoms, 
either bosonic or fermionic, is a point already consid¬ 
ered [19, 20]. In particular the possibility of constraining 
the energy-momentum relation resorting to cold atoms 
has already shown us that this kind of systems could 
open up new landscapes in the context of gravitational 
physics [21]. Of course, this last topic in the context 
of phenomenology of quantum gravity leads us to ask if 
there are additional low-temperature effects that could 
be used as trackers for new effects. This question implies 
the quest for the detection of these kind of effects in the 
realm of condensed matter physics, i.e., a broadening of 
the current attempts. Clearly, another low-temperature 
effect is the phenomenon of superfluidity [22] and, there¬ 
fore, we wonder if this case could offer a new window for 
our search. In the present work we explore this situation, 
namely, the possibilities that superfluidity has to offer in 
the context of quantum gravity phenomenology. 

Having stated our goal we must discuss, though briefly, 
the physics behind the emergence of viscosity in a flow. 
Concerning the phenomenon of superfluidity the first ex¬ 
perimental results can be found in the work of Kamerling 
Onnes of 1911 in which he detected that if cooled below 
2.2 “AT He did not contract but rather expand [23]. The 
current work has been able to provide a coherent pic¬ 
ture to the subjacent Physics [24, 25]. Bose-Einstein 
Condensation (BEC) is also connected to the presence 
of very low temperatures. Fritz London [26] put forward 
the idea of a connection between these two effects assert¬ 
ing that the transition from He I (the high temperature 
phase of liquid helium) and Hell (the low temperature 
phase) should be considered an example of a BEC. Tak¬ 
ing into account London’s idea and, joining it to the pre¬ 
vious work BEC-quantum gravity phenomenology, once 
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again, we are confronted with the question about a pos¬ 
sible use of superfluidity as a tool in our quest. 

The concept of elementary excitations in the realm of 
superfluidity was first introduced by Landau [27] (within 
the two-fluid model proposed by Tisza [28]) as a core fea¬ 
ture in the description of the behavior of He II. Landau 
asserted [29] that the normal fluid (the non-superfluid 
component) could be regarded as a dilute gas whose com¬ 
ponents are weakly-interacting elementary excitations 
which move in a background defined by the superfluid 
component. Along these ideas the phenomenon of super¬ 
fluidity appears when the velocity of the corresponding 
flow lies below a certain threshold value given by 



In this last expression e(p) denotes the energy of an 
elementary excitation and p the corresponding momen¬ 
tum. 

In the experimental realm the quest for this critical 
velocity has been carried out in a sodium-BEC, and the 
results show a possible velocity threshold located around 
the value of 1.6 mm/s [30]. Sodium is a system that can 
be condensed and the speed of sound in it has already 
been detected [31, 32]. 

At this point we may now state clearly the ideas con¬ 
tained in the present work and why they have been cho¬ 
sen. The main purpose is to obtain a prediction for the 
critical velocity for a BEC in which the relation energy- 
momentum of the particles of the gas has been deformed 
along some possibilities contained in several quantum 
gravity models. The system is a sodium gas trapped 
by an isotropic harmonic oscillator potential. The rea¬ 
son for this particular choice stems from the fact that 
several of its properties have already been detected and 
measured, for instance, evidence for a critical velocity [30] 
and speed of sound [31, 32]. In other words, our proposal 
is the following one: use a BEC-condensed sodium gas 
and measure the region in which superfluidity is present. 
Compare the size of this region against the theoretical 
prediction here obtained and deduce the corresponding 
bound for our parameter containing the breakdown of 
Lorentz symmetry. Up to now, superfluidity has not 
been considered a relevant element in quantum gravity 
phenomenology and the present work shall be considered 
as an analysis of the perspectives that this topic could 
offer. 


II. MATHEMATICAL MODEL 

Even the simplest mathematical model for a BEC (the 
Gross-Pitaevski equation) trapped by a harmonic oscilla¬ 
tor potential has no analytical solution, yet [33]. This last 
fact implies that in this topic we must resort to approxi¬ 
mation methods, etc. In our case, the theoretical model 


will describe the features of the BEC resorting to an ap¬ 
proximation method in which the presence of interactions 
among the particles produces a change in the frequency 
of the trap rendering a smaller value than the one pro¬ 
vided by the trap, for the case of repulsive interactions. 
This assumption will allows us to calculate the energy of 
the ground state and of the thermal cloud. The energy 
and momentum of the elementary excitations, according 
to Bogoliubov ideas, are a function of the energy of the 
excited particles, and this parameter can be computed 
from our assumptions. Finally, the critical velocity is 
deduced as a function of our phenomenological variables 
and compared against the measurement readouts. From 
this comparison an upper bound for our model values will 
appear. 

From a fundamental point of view our mathematical 
model can be defined by an A^-particle Hamiltonian the 
one in the formalism of second quantization is [34] 


H = df 


Uo 


2m 

+V{r)ip^{r,t)tp{r,t) 


—(f, (f, t)V>(r, t)V;(r, t) 
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In this Hamiltonian and tp{r,t) represent 

bosonic creation and annihilation operators, respectively. 
It is valid only at low energies and momenta and implies 
that the interaction among the particles is, as usual, cod¬ 
ified by the scattering length parameter a, i.e., Uq = 
. The corresponding trapping potential (U(r)) is 
an isotropic harmonic oscillator whose frequency reads 
oj. Moreover, there are N particles in the gas, each of 
them with mass m, the volume occupied by the system 
is V. 

The mathematical assumptions in the present model 
read: 

(i) Only two states are populated, namely, ground 
and the first excited state. We may provide to this as¬ 
sumption a realistic physical meaning recalling that for 
a bosonic system, with chemical potential p and energy 
levels of single-particle e, the occupation number in ther¬ 
mal equilibrium is given by [35] {j3 = 1 /{kT)) 


< >- g(e-M)/3 _ 1 ■ ( 4 ) 

Clearly, it is a monotonic decreasing function of e, and 
this feature justifies the present assumption. 

(ii) The Hartree approximation will be employed for 
the mathematical description of the two occupied states. 
In other words, the ground state of the interacting sys¬ 
tem is deduced by a Ginzburg-Pitaevski-Gross energy 
functional [36], and it entails that the ground state wave- 
function corresponds to the case of a harmonic oscillator 
situation but the frequency is modified due to the fact 



3 


that the system has a non-vanishing scattering length 
[37], such that the fundamental length parameter reads. 


R = 




(5) 


In this last expression I is the radius related to the trap 
given by the isotropic harmonic oscillator of the trap 


I = 



( 6 ) 


Of course, we end up with an effective frequency 


h 

mR"^' 


(7) 


The experimental conditions entail R > I [30] and, 
therefore, oj < to. 

The order parameter of those particles in the ground 
state is 


/ ^(0) 

r J.2 




( 8 ) 


Here fV(o) denotes the number of particles in the lowest 
energy state. The presence of a non-vanishing scatter¬ 
ing length implies that in the ground state not all the 
particles can have zero-momentum, the reason for this 
lies in the fact that the two-body interaction mixes in 
components with atoms in other states [35] and 


= (12) 

From the symmetry of the system we may conclude 
that an excited particle can have vanishing momentum 
along the x and y axes but a non-zero one in the z- 
direction, or vanishing momentum along the z and y axes 
but non-zero one in the x-direction, or, finally, vanish¬ 
ing momentum along the x and z axes and larger than 
zero along y-direction. The symmetry of our trap and 
of the scattering length entail that one third of the ex¬ 
cited particles will have non-vanishing momentum along 
the x-axis, one third along the y-axis, and the remaining 
third along the z-direction. Mathematically this corre¬ 
sponds to the following expressions 



8 In iNa^x^^ r n 


(13) 


Here x^^'> = x, x^^'> = y, and x^^'> = z. 

Of course, (13) must be related to the total number of 

particles in excited states = ^N^ T^); ^ condition 
that becomes [35] 


N 


(e) = 




d^r. 


(14) 


The three expressions in (13) will be used for the defi¬ 
nition of the wavefunction (4'(i)(r)) of the thermal cloud, 
i.e., 


Clearly, 


7V(o) = N 


, 8 

“ ttH . ■ 


^(i)(7') =V'(i)(7') + V’(ij(rO+V'(ij(r). (15) 
(9) 

III. ELEMENTARY EXCITATIONS AND 
DEFORMED DISPERSION RELATIONS 


N(o) = J ( 10 ) 

V=—R^. (11) 

Concerning the thermal cloud, the core of this part is 
also comprised by the Hartree approximation. We as¬ 
sume that all excited particles are in the same state and, 
since the temperature is very low, it corresponds to the 
first excited state of a particle trapped by a harmonic 
oscillator with a frequency given by (7). If 4'(i) denotes 
the wave function of the thermal cloud, the wave- 
function of the first excited state of a single particle in 
our effective trap, and fV(g) the number of particles in the 
cloud then 


Having stated our assumptions we proceed to compute 
the critical velocity [22]. The deduction of the energy 
of an elementary excitation and of its corresponding mo¬ 
mentum requires the knowledge of the energy of a single¬ 
particle in the first excited state [34]. The thermal cloud 
contains particles in the first excited state of an isotropic 
harmonic oscillator whose frequency is (7) therefore the 
energy of an excited particle is given by this assump¬ 
tion and easily calculated as a function of the effective 
frequency of our variational procedure 

g(o) = ^hib. (16) 

This is the case in which no deformed dispersion re¬ 
lation has been considered. Introducing this quantum 
gravity parameter we have 
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e = 


^ ~ n 

— nu + ap 


(17) 


a=( 


47r\ 2/3 r 


.3-) 






47r\ 1/3 in 


a. 


(23) 


In this last expresion a and n are parameters stemming 
from the considered quantum gravity model [18]. 

According to the ideas of Bogoliubov [34, 38] the en¬ 
ergy of an elementary excitation, here denoted by e, is 
a function of the energy of the excited particles of the 
BEC, namely, 


e = 



2iV[/(o) . 


(18) 


The contribution to the energy of all the elementary 
excitations turns out to be [34, 38] 


According to Landau [29], in order to find the criti¬ 
cal velocity we must now deduce the momentum of this 
elementary excitation. These physical variables, which 
define the normal component of the fluid fluid, can be 
regarded as a bosonic gas whose components are weakly- 
interacting and moving in a region in which a constant 
potential exists, and this potential is defined by a mean 
field approach [34]. According to this interpretation we 
may rewrite (22) in the same form as in the case in which 
our BEC is a homogeneous one [34]. In other words, take 
(22) impose the condition a = 0 and compare the result 
against 




2NUi 


( 0 ) 


V 


e < he > . 


(19) 


€ 


2m V V 


(24) 


We have defined < he > as the occupation number of 
the elementary excitations with energy e. The relation 
between the occupation numbers of particles and elemen¬ 
tary excitations is [34] 


This last argument allows us to deduce the wavenum¬ 
ber related to our elementary excitation and, therefore, 
its momentum. 


< he >= 


< rie > 
1-1- < rie > 


( 20 ) 


Since our model has to be consistent with the present 
experimental technology at this point we resort to the 
laboratory values related to the detection of a critical 
velocity in a sodium condensed gas [30] in which the 
occupation number of the particles in the first excited 
state fulfills the condition Af(g) ~ 10^ > 1. Therefore, 
< >= 1. In addition, < ne^^y >= 0, Vi > 1. Indeed, 

we have considered that the thermal cloud is comprised 
by particles which occupy only the first excited state, in 
other words, < >= 0, Vi > 1. Introducing this 

condition into (20) leads us to the aforementioned result 
for the occupation number of the elementary excitations. 

The order of magnitude of this deformed dispersion 
relation has to be very small, otherwise, it would have 
already been detected. This means that we must expect 
the fulfillment of 



Resorting to Landau criterion (2) we obtain that the 
critical velocity is given by 


1 h /25 /47r\2/3 207rNa 
1 /47r\i/3 

^TIoVYJ 



(27) 


IV. CRITICAL VELOCITY 


>> ap”. (21) 

We now cast (18) in a different form, and for this we 
resort to the effective volume V = use (5), (6), 

and (7) and keep only terms linear in a. The final result 
reads 


e 


/47r\i/3 /25 /47r\2/3 207rVa 

vYJ mE2/3 V TVYJ Vi/3 


a I AttN a 1 

~^2y 51/1/3 y«/3> 


( 22 ) 


At this point we proceed to check our model, and in 
order to do this we consider the case in which a = 0. The 
experimental parameters [30] to be used are: (i) a critical 
speed of = 1-6 mm/s] (ii) the number of particles 

in this experiment has a minimum of V = 3 x 10® and a 
maximum of V = 12 x 10®, and for the evaluation of our 
expression we will take the arithmetic average, i.e., N = 
7.5 X 10®; (iii) the effective volume is that of an ellipsoid 
whose axes are = 45 x 10“®m and li = 150 x 10“®m 
such that V = and, finally, (iv) a scattering length 

a = 2.75 X 10-®TO. 
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Introducing these values into (27) (setting a = 0) im¬ 
plies 

y(m) _ 2^ gg rnm!s. (28) 

The reported critical speed is [30] 

= 1.6 mrajs (29) 

The ensuing error is less that 18 percent 

= 0.179. (30) 


a choice has to be made in connection with the number 
of particles since this physical variable changed from ex¬ 
periment to experiment; for the sake of concreteness we 
consider the highest value, i.e., fV = 12 x 10®, taking the 
lowest case {N = 3 x 10®) does not modify the order of 
magnitude of the ensuing bound. 

The measurement readouts [30] are given up to units of 
tenths of mm/s and, in consequence, the smallest scale 
can be considered as an approximation for the experi¬ 
mental error [39]. In other words, if Av denotes the ex¬ 
perimental error of the measuring device, then the afore¬ 
mentioned argument implies An ~ 0.1 mm/s. The ex¬ 
perimental error has to be equal or larger than the term 
containing the effects of the breakdown of Lorentz sym¬ 
metry, this phrase means for this situation 


In other words, our model provides a very good de¬ 
scription of the experiment and, hence, the analysis of 
the deformed case within the present framework seems a 
reasonable assumption. 


V. CONCLUSIONS 

In the present work the analysis of the options that su¬ 
perfluidity could offer in the context of quantum gravity 
phenomenology has been done. The model has been a 
Bose-condensed sodium gas trapped by an isotropic har¬ 
monic oscillator in which the energy-momentum relation 
for the particles has been deformed along the propos¬ 
als emerging from some quantum gravity models. Af¬ 
terwards, along a perturbative approach, the energy and 
momentum of the elementary excitations generated by 
the particles in the thermal cloud have been calculated. 
Finally, we introduce these last two physical parameters 
in Landau criterion associated to superfluidity and hnd 
that the breakdown of Lorentz symmetry, in the form 
of a deformed dispersion relation, implies a modification 
of the region in which, for a sodium BEC, superfluid¬ 
ity may exist. If a > 0, then the aforementioned region 
grows (compared to the case in which a = 0), whereas, 
if a < 0, then this region becomes smaller. Indeed, the 
allowed superfluidity velocities are those falling into the 
interval (0,V(crit)). and, clearly, this interval becomes 
larger for a > 0. 

In relation with a bound for our phenomenological pa¬ 
rameter a, a fleeting glimpse at (27) tells us that we must 
first choose a value for n. As an example we take n = 1 (a 
condition that implies that a has units of speed) and con¬ 
sider the experimental values related to the evidence of a 
critical velocity in a sodium condensed gas [30]. Clearly 


14 X 10 ® mm/s > a. (31) 

Additional cases (n = 2 does not imply new physics, 
it entails only a redefinition of the inertial mass) can 
be analyzed in the same manner. The present argument 
tells us that we may deduce an upper bound for deformed 
dispersion relations associated to the structure given by 
(1), and this expression reads 


Av > 


47r\ 1/3 


^(i) 



Let us comment that in the present model the de¬ 
formed dispersion relation has been introduced only in 
the context of particles (see (17)) but not in connec¬ 
tion with the corresponding elementary excitations also 
called quasi-particles, expression (22)). Clearly, an ad¬ 
ditional possibility is the introduction of the breakdown 
of Lorentz symmetry at the level of the kinematics of 
the quasi-particles. This second option implies the in¬ 
troduction of a second pair of phenomenological param¬ 
eters, since the quantum gravity modihcations could be 
particle-dependent. This case can be, without any fur¬ 
ther problem, be considered in the present framework. 

Summing up, we have shown that superflnidity in 
sodium-condensed gases offers a new window in the 
realm of quantum gravity phenomenology and that the 
present experimental results are enough to deduce some 
rough bound for the involved parameters. 
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